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Abstract 

We consider moving particles in media with nonlinear friction and drive them by an 
asymmetric dichotomic Markov process. Due to different energy dissipations, during 
the forward and backward stroke, we obtain a mean non-vanishing directed flow of 
the particles. Starting with the stationary velocity distribution, we calculate the 
stationary current of particles, the variance and the relative variance in dependence 
on the degree of nonlinearity of the friction, on the asymmetry and for different 
strengths of friction. In two dimensions the particle performs diffusional motion, if 
in addition the direction of the asymmetric drive changes stochastically. 
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1 Introduction 



The occurrence of directed flows of matter or single particles, in case of mean 
vanishing forces, is always connected with a break of symmetry in the consid- 
ered medium or device. Since the pioneering work of H. Purcell [1] a lot of work 
has been devoted to this problem. Purcell found that simple devices require 
more than one degree of freedom by which they can temporally store energy 
sequentially in an asymmetric way. Later on, theoretical work on ratchets 
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has elaborated the many different combinations of possible constructions be- 
tween nonlinear periodic potentials, noise and forcing to obtain non- vanishing 
flows of particles if the averaged sum of forcing vanishes [2"I3THIoTST7] . It was 
P. Hanggi who put forward the notion of Brownian motors [5115115] to point 
out the possibility to find engines in the presence of a strongly fluctuating 
environment. Most strikingly, in many cases this noisy environment can even 
act as the source of energy which drives the engine (the interested reader is 
referred to [10] for Peter Hanggi's extended list of publications on Brownian 
motors and the various different topics he has investigated). 



In this paper we study particles moving in a medium with nonlinear friction 
[TT|12||13|[Ti] . Therefore, the rate of dissipation will depend on their velocity. 
Achieving different values of velocity will result in a different speed of dissipa- 
tion. We will drive such particles by random forces with a vanishing average. 
We assume a two state random telegraph signal or a dichotomous Markov 
process as force acting on the particle [T5lll6|ll7|[T8] . Crucially, we will assume 
that this force acts asymmetrically, i.e. its two values in the two directions 
will be different. Such an applied asymmetric force elongates the velocity of 
the particle differently during the forward and backward stroke. Hence, it will 
dissipate the transmitted energy during both strikes differently. For this rea- 
son the particle is able to perform a directed motion even in the case if the 
applied drive vanishes on average. 



We will assume that the random drive is a time homogeneous stochastic pro- 
cess. Hence, the problem becomes stationary in the long time limit and we 
can use tools from the theory of stationary stochastic processes. As a result 
we calculate the mean velocities and the variance of the velocities for the one 
dimensional case. In two dimensions, we add a stochastic rotation of the di- 
rection of the applied asymmetric forcing and show first numeric results for 
the diffusional motion of the particle. 



We mention that similar problems have been discussed within the framework 
of vibrational dynamics by Blekhman (see e.g. [19]), who has used an alterna- 
tive approach in the case of a fast periodic driving to find asymptotic velocities 
of gliding particles in many different situations. This approach is based on the 
separation of time scales and averaging over the fast time scale of the driving 
force. Similar techniques have also been used as approximate solutions to elab- 
orate instabilities in driven oscillators, see e.g. the Kapiza- (Magnus) pendulum 
in [20]. An application to our stochastic driving is also possible, but here we 
rely on techniques using the stationary probability density for the particle's 
velocity. 
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2 Nonlinear friction and asymmetric forcing 



We study singular particles with unit mass m = 1 which move in a fluid. 
Therein the particles experience a frictional force f(v) which damps the mo- 
tion. To compensate the damping, the particles are forced by additive tem- 
poral forces £(t). We assume that this force possesses a vanishing mean and 
stationary correlations 

<£(*)> = o, mt(t)) = K(t-i?) . (i) 



Thus the velocity v(t) of the particles obeys the following dynamical equation 

v ( t ) = f (v( t )) + at) . (2) 

The time-dependent force £(t) drives the system permanently out of (thermal) 
equilibrium. The velocity distribution is non-Maxwellian. We underline that 
there is no spatial dependence in the description ([2]). This typical ingredient of 
Brownian ratchets, namely spatially periodic gradients, is absent in our model. 
Thus the rectification process of ([2]) differs from the usual Brownian ratchet 
mechanism [3J. As outlined, symmetry is broken due to the asymmetric forcing 
realizing an asymmetric energy dissipation by the nonlinear friction. 

As driving force we will consider a random telegraph signal £(£) which is 
a time- homogeneous Markov process. It can take two values with 
constant transition rates between the two states. k + denotes the rate 

of passage from A + to A_ and analogously, k_ denotes the rate of passage 
from A_ to A + . The resulting time correlation function reads 

(£(*)£(«))= {A+ (~ A Tk» K ~ ■ e- (K++K - )lt - sl ■ (3) 

[ K i -+- K_ ] 



Particular interest will be paid to unbiased dichotomous Markov processes, 
i.e. processes with a vanishing mean (£(£)) = 0, which leads to 

A + k_ + A_k + = . (4) 



The condition reduces the amount of independent parameters to three. We 
therefore introduce the following set of three independent parameters of the 
unbiased DMN process 

A=\A-\ + \A+\, t = t_+t + , 



< p = 



min(|A_| 


• 1 


A+\) 


max(L4_ 


1 1 


\A + \) 
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where the mean waiting times t± for the states A± relate to the rates as 
t± = l/k±. The parameter A measures the strength of the DMN process, r, 
the mean time of one cycle, is the characteristic time scale of the DMN. It 
differs from the correlation time tq = (k + + The parameter p will be 

the desired parameter controlling the asymmetry of the driving. Symmetric 
driving is located exactly at p — 1. High asymmetries, i.e. when the stroke 
amplitudes and \A + \ differ much, let p tend to zero. 

We can express the variance of the DMN in terms of the new parameters, 
resulting in 

varm = (e(t)} = ^- 2 . (6) 

The variance increases with the driving amplitude A. This is not rather sur- 
prising. However, while looking at the variance profile for a varying asymmetry 
p, one observes, that the variance reaches its maximum when the driving is 
totally symmetric. On the other hand, growing asymmetry (i.e. lower values of 
p) reduce the variance. In the asymmetric limit p — > 0, the variance vanishes 
completely. 

In particular we are interested in the effect of nonlinear friction. Therefore it 
will be convenient to introduce a model of friction that can model nonlinear- 
ity via a parameter, in order to easily control the impact of nonlinearity on 
the system (T5]). Prominent types of friction, are the so called Stokes friction 
(f(v) = —j s v), which is typical for small particles moving at low velocities 
through viscous fluids, and the "quadratic drag force" (f(v) = — jq\v \v ), corre- 
sponding to situations of objects moving at relatively large velocities through 
fluids, a widely occurring situation in aerodynamical engineering. Another 
prominent type of friction is the so called Coulomb friction (f(v) = — 7 c o0- 
This form of friction characterizes sliding frictions, the form of friction that 
occurs, when two surfaces slide against each other. This situation is sometimes 
called "dry friction" . The respective friction coefficients 7$ depend specifically 
on the characteristics of the fluid and of the corresponding particles. 

Moreover, Stokes friction is a linear friction model, while the other two models, 
quadratic and Coulomb friction, are nonlinear. Our model of friction f(v) 
should be able to model more general forms of friction, but at the same time 
able to reproduce the mentioned three prominent models of friction as specific 
limit cases. An immediate candidate for such a model is apparently 

f(v) = - 7 H n sign(t;) 7 > 0, n G R+ . (7) 

This friction model introduces two parameters, namely the friction coefficient 
7 and the friction exponent n, allowing to model a wide range of relevant types 
of friction beyond the introduced cases (n = 0,1,2). The friction of general 
fluids, that are omnipresent in small scaled biological systems, may therefore 
be properly modeled by (J7]) (see e.g. [T9]). 
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In the other limit, for n — > oo, (J7J) reduces to 

!+oo for v > 1, 
for |M| < 1, (8) 
— oo for v < —1. 

The motion of the particle is constraint to the velocities between —1 and 
1, since the damping friction forces flH]) approaches infinity for \v\ > 1. In 
the region \v\ < 1 the particle moves freely, without damping, while at the 
height of the barrier sites \v\ = 1, the particle gets reflected back to lower 
velocities. This is in some sense similar to a particle moving in a potential 
well, with infinite slopes of the respective walls. Note that in this case there 
is no dependence on any friction constant 7 anymore. 



3 Stationary currents: Analytically tractable models 



In general, it is impossible to obtain an exact solution of, for instance, the 
stationary density of a system, if an arbitrary form of noise, especially colored 
noise, is considered. However, the simple structure of the DMN, as a Markovian 
two-state process, enables us to obtain exact and explicit results in terms of 
stationary distributions [TD r )18f2ip^] . The stationary solution P st for Eq.(J2]) 
reads 



P st (y) = N\i + (v)i_(v)\ 



x exp 



where the functions i±(v) are defined as 



k + I i + (s)ds — k_ / i_(s)ds 



(9) 



i±lv) 



-7|v| n sign(t>) + A± 



(10) 



As one can see with Eq.()2]) and (J7j), the support of the stationary distribution 

is a compact interval with v + = ^A + /j and V- = — ^A~ZJj, and 

beyond this interval the stationary distribution has to vanishes. 



In order to compute the integrals in Eq.(9), we look at the general expression 
for arbitrary a, C,x,n G R. Notice that we can always write 



ax n +C ' 



1 



ax n + C C 



E 

i=0 



, ax . ^ 



as far as < 1 holds. 
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Integrals of this expression can be expressed by the the standard classical (or 
Gaussian) hypergeometric series 2 Fi. In its general form, it is defined as 



2F1 [a, b, c,z} = J2 



(a) n (b) n z r - 



n=0 vCjn U\ 



for a,b,c, z G C 



(12) 



Integration of Eq.(ll) yields 



r 



dx 



= ^ 2 Fx 



ax n + C C 



-,1,1 + -,-*" 

n n C . 



(13) 



and consequently 



/-u 
i±(u)d 



f = < 



4* 



2F1 



± 



1 



,1.1 + 
n n A± 



v>0, 



v < 0. 



Stationary moments can be computed straightforwardly. However, in general it 
will be difficult to give an analytic closed formula for the stationary moments. 
That is why we will first look at special cases and discuss the general case of 
arbitrary n G M + in the context of simulations later on. 

Let us start with the dry friction limit (n = 0), i.e. 



v = -7 sign(v) + £(t) 



(14) 



Notice that for the case 7 > max(\A_\, \A + \) no motion can occur at all, since 
the critical force in order to induce motion is ±7. In this case the stationary 
distribution P s ° takes the form 



(15) 



where the 5 function denotes the Dirac distribution. 



To exclude this situation, we require 7 < min(\A_\, \A + \). Then the critical 
friction force is overcome in both situations, in the state A + as well as in the 
state A_ of the DMN 
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The functions i±(v) simply become constants which results in 



-Ni 



exp 


K + V 










:-7 + a+) ( 


-1 + A-) 



v > 0, 



(16) 



exp 



K + V 



K_V 



7 + A+ 7 + A_ 



( 7 + A+) (7 + A_) 



v < 0. 



The stationary distribution P s ° t is an exponential function for every half line 
{y > 0, v < 0). Note that for (fl6|) to be normalizable, the exponentials should 
converge to zero in the limits of infinity v — > ±oo. Therefore it is sufficient 
to show, that the exponentials on every half line are decaying, which is in 
accordance with our assumptions A + > 0, A_ < 0, re±, 7 > and 7 < 
min(|A_|,|A + |). 

If the driving is asymmetric, we get two different limits for the distribution 
function at the origin v — 0, since 



lim P» 

D->0+ 



N 



[—j + A+) (-7 + A. 



(17) 



and the left side limit yields 



lim P» 



A' 



(7 + A + )( 7 + A_) 



Thus, the stationary distribution is discontinuous at v = in the dry friction 
limit (n = 0). 



For the normalization constant iV we get 



iV = — — 



7K 
2" 



(19) 



with k = k + + k_. The first moment, as the mean average velocity and sta- 
tionary current J = (v), gives 



(v) 



A- +A, 



K 



(20) 



Hence, in the dry friction limit, the particle moves in the direction of the 
stronger of both strokes A±. If \A + \ > |A_|, the particle moves into the positive 
direction while in the opposite case it moves backwards. Remarkably, (v) does 
not explicitly depend on the (critical) friction force 7. No matter how strong 
the stiff friction force 7 is, as long as 7 < min(\A_\, A + \) is ensured, the 
average velocity (v) does not depend on the friction. Expressing the mean 
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A= 0.5 
A= 1.0 
A= 3.0 
A= 5.0 
A=10.0 




Fig. 1. Case of dry friction: current versus asymmetry parameter p for different 
stroke amplitudes A. The current shows a maximum for optimal asymmetry 



velocity (120]) in terms of the three effective parameters of the driving p, A and 
r, we obtain 

_p(p ~ 1) 



( v ) = At- 



(21) 



(p+l) 3 ' 

where we assumed, without loss of generality, |A_| > \A + \. Remember that we 
have per definition < p < 1. Thus, in this case the current is negative. For a 
symmetric driving (p = 1) the velocity vanishes exactly. However, the velocity 
vanishes also if p — > 0. That is why too much asymmetry is not a successful 
way to obtain the optimum speed in the dry friction limit. Therefore we expect 
a maximum of the mean velocity within the region p G [0, 1]. This result is 
confirmed by the numerical simulations shown in Fig.(l). 

We now turn shortly to the case of Stokes friction (n — 1 or f(v) = —jv) and 
will take a look at the mean velocity. For the stationary distribution P} t we 
compute 

Pl t (v) = N \-yu + A + fr~ l |-7u + A-fr' 1 . (22) 
To determine the first moment J = (v), we get 



(v) 



~1~ ' 1 1 

v (— 7f + A + )~ (— 7t> + AJ)~ dv 



r[? + 1] 



{A+K. 







(23) 
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Thus we have obtained, that in the linear case, there can never be a current 
for vanishing mean forcing. The ratchet mechanism does not work for linear 
friction. 

In the limit of n — > oo, there is no current as well. Since the relation (8) 
holds, there is no friction forces in the interval v G [—1,1]. Here P™(v) is 
constant. Hence, the first moment vanishes and for the second moment we get 
(v 2 ) = 1/3. 



4 Stationary current: The general case 

In general cases (n G the analytical treatment is difficult. Therefore we 
will turn to numerical measures and our primary focus will lie on the mean 
velocity (v) which equals the stationary current. In addition, we will look at 
the variance var {v} and at the inverse relative variance 

coh{v} = 2 [ V) \ (24) 

which might serve as a measure of coherence of the transport [23]. High co- 
herency means, that one has a low spread in the transport process compared 
to the unidirectional drift. 

The impact of the nonlinearity parameter n on the mean velocity (v) is dis- 
played in Fig. ([2]) (with n G [0,4]) for different realizations of the friction 
coefficient 7 and the driving asymmetry p. In all these cases we observe a sin- 
gle current reversal exactly at the linear friction regime n = 1. In other words, 
there is no transport for linear friction. This is exactly the result we have found 
analytically. In particular, sublinear {n < 1) and superlinear (n > 1) friction 
results in opposing velocity directions. While for sublinear friction, one ob- 
serves a negative current, we get a positive current for the superlinear friction 
regime. Note that, without loss of generality, we have assumed |A_| > \A + \. In 
the opposite case (\A + \ > |A_|), one gets the same picture with inverse sign. 
While the sublinear regime (dry friction regime) results in the direction of the 
stronger amplitude, the superlinear regime favors the opposite direction. In 
the dry friction limit, we observe a different behavior. If 7 < max(|A_|, \A + \) 
holds, the current becomes a non-zero value. Otherwise the current vanishes. 
This is exactly what we predicted already in the analytical study of the dry 
friction limit n = 0. In the superlinear case the stationary current does weakly 
depend on the friction. In the limit of large n, the current will disappear. 

The variance of the velocity is shown in Fig. ([3]). When the critical force cannot 
be overcome by the driving, the variance again vanishes in the dry friction 
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n 

Fig. 2. Current J versus friction exponent n for different friction coefficients 7 (top) 
and driving asymmetries p (bottom) |29j . 

limit. In this situation the variance increases with growing friction exponent 
n. In the case, where the critical force can be overcome by the driving forces, 
we observe a finite variance for n — 0, which decreases with increasing friction 
exponent. 

The coherence is presented in Fig. ([3]). Since the current J vanishes for linear 
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Fig. 3. Variance of velocity for different friction coefficients 7 (top) and coherence 
(Eq. (|24p ) for different driving asymmetries p (bottom) versus friction exponent n 

friction and at the same time the variance keeps finite, we observe that the 
transport coherence tends to vanish. We also see, that for superlinear friction 
(n > 1), the dependency of the coherence on the exponent n is very weak. 
coh{v} almost does not change with varying friction exponent n. Whereas 
we observe a rather strong dependence on the exponent n in the sublinear 
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regime. In the limit n — > 0, the coherence increases exponentially. Thus, the 
dry friction regime shows much better coherence than the friction regime, 
which occurs in fluids. 




Let us now investigate the asymmetry parameter p. In Fig. (j4]) one observes in 
all cases a single maxima at symmetry percentages below 0.5, varying slightly 
for different parameter settings. In accordance with our previous considera- 
tions, the current vanishes for a symmetric driving as well as in the strongly 
asymmetric case. This is elucidated by the specific values of driving amplitudes 
A± and waiting times r± in this limit 

limA_ = —A, UmA + = 0, (25) 

lim r_ = 0, lim r + = r . (26) 

p^O p->0 v ' 

The average waiting time r_ for the state A_ vanishes. Consequently, no nega- 
tive force A_ acts. On the other hand, no positive force acts as well, according 
to (25). Hence, no forces at all are acting, if p = 0. 

The same line of argument shows, that the variance also vanishes in the asym- 
metric limit. Indeed, the results confirm our ideas (Fig. (5)). Thus, starting at 
zero for p = 0, the variance, in any case, increases, until it reaches its maxi- 
mum for p = 1. This can be understood in terms of the driving £(£), which has 
already shown this behavior (Eq. (6)). It seems reasonable that the variance 
of the state variable v follows qualitatively the variance of its driving. 
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Fig. 5. Variance (top) and coherence (Eq.(24);bottom) for different friction coeffi- 
cients 7 versus driving asymmetry p [29\ . 



The behavior of the coherence (Fig.(jS])) shows similarities to the behavior of 
the current J. One observes a single maximum, which this time is shifted to 
higher asymmetries, i.e. to lower values of p, compared to the location of the 
maxima of the current. For the regarded cases, the peaks of the current are 
in regions where p > 0.2 holds, while the peaks of the coherence are located 
below that barrier. We observe a much better localization of the maxima of 
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the coherence, compared to those of the current. 



5 Diffusion in the two dimensional case 



Eventually, we present results for the two dimensional case. We assume that 
the asymmetric drive £ (t) acts on the particle along a certain variable direction 
in the two dimensional space. This direction is given by the unit vector with 
Cartesian coordinates 

■ fcSS) <*> 

and 0(i) is the present angle between the direction of the drive and the in- 
coordinate. Along this axis, the particle changes v(t) with positive and nega- 
tive values according to Eq.(l) and 

^-r(t)=i?(t) while v(t) =v(t) ■ e v (t) . (28) 

Additionally, a Gaussian white noise source, labeled rj(t), rotates the moving- 
axis of the particle. As equation for the angle Q(t), we formulate 

£e(f) = - V (t), (29) 
at v 

which assumes that the increment of the rotation angle scales with 1/v. The 
noise rj(t) is defined by (r](t)) = and (r}(t)v(t')) = 2D v 5(t - £'). 

For the following simulations, we fixed the friction exponent n = 3 with \A + \ > 
\A_\, the intensity of the noise D v = 0.1 and all other parameters as stated in 
EDI. 



The motion of the particle shows a diffusive behavior after a crossover. Due 
to the noisy rotations of the asymmetric force, no preferred direction exists. 
Typical trajectories are presented in the inset of Fig. OH]). Simulations show 
that the velocity distribution densities have, for not too small r, a vanishing 
derivative at v — >■ and possess two peaks according to the asymmetric drive. 
This can be expected in regard to the one dimensional case ([23] and Eq.(9)). 
In the two dimensional velocity space, two circles with high probabilities occur. 

For large time scales, we assume ((r(t) — r(t )) 2 ) = 4tD e ff. The resulting be- 
havior of the diffusion coefficient (Fig. (JB])) can be explained by the trajectories 
in the inset. They illustrate the impact of the small force A_ which acts most 
of the time for p = 0.1 and which leads to a stronger influence of the angular 
noise, compared to the symmetric case p = 1. 
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Fig. 6. Diffusion coefficient of the two dimensional motion versus driving asymme- 
try p for several and different friction exponents n; spatial trajectories for two 
different driving asymmetries p in the inset [29] . 

As a possible explanation, we use the Taylor-Kubo-relation [251126] 




In simulations we have inspected that the velocity correlations CW(r) increases 
for stronger symmetries p in the driving due to an increased standard devia- 
tion of the velocities. Therefore the diffusion coefficient -D e // is maximal for 
symmetric noise. 

The slowly growing diffusion for small friction exponents and small driving 
symmetries p (Fig.(jBj)), can be explained by the nearly vanishing velocity u_, 
which mostly acts in this limit and depends on A_. We also observe that larger 
noises D n lead to smaller diffusion coefficients. This behavior has already been 
studied in similar systems [25|27] and is caused by a straight movement for 
small angular noise. The analytical treatment of the diffusion coefficient, for 
the one dimensional DMN process, has been recently studied for a symmetric 
driving [28] and corresponds to our results for p — 1 in the limit of small noises 



6 Conclusion and outlook 

We have shown that directed motion in nonlinear media is possible, if particles 
are driven by an asymmetric force. The force can have external or internal 





(30) 
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origin. Small additive noise would only partially destroy the obtained effect as 
was seen in calculation of the stationary distribution densities [18] . In the one 
dimensional case, it results in a smoothing of the densities but the averaged 
motion will survive. Diffusion will be obtained with properties similar to the 
recently studied case with symmetric drive [28] . 

While the analytically and numerically studied one dimensional case is well 
understood, the behavior of particles in two dimensions remains still a chal- 
lenge. Also further energetic consideration will be helpful in order to estimate 
the efficiency of the proposed ratchet. Future fields of interest will include also 
alternative two dimensional models for DMN processes and the influence of 
coupling forces to external obstacles or mutually to similar objects. We also 
see a great challenge in the consideration of different friction terms or peri- 
odically changing friction forces in order to simulate rhythmic movements in 
biological systems. 
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